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ABSTRACT 

We present a comprehensive spectral analysis of black hole X-ray binaries, LMC X-1 
and LMC X-3, based on BeppoSAX observations. We test both the multi-color disk 
plus power law (MCD-I-PL) model and a newly-developed Monte-Carlo simulation- 
based model for a Comptonized MCD (CMCD) with either a spherical or a slab-like 
corona, by comparing the inferred parameters with independent direct measurements. 
While all models give an adequate description of the spectra, we find a significant 
discrepancy between the MCD-t-PL inferred X-ray-absorbing gas column density and 
the absorption-edge measurement based on dispersed X-ray spectra. The MCD-I-PL 
fits to the LMC X-1 spectra also require a change in the inner disk radius during the 
BeppoSAX observation, which may be due to the nonphysical effects inherited in the 
model. In contrast, the CMCD model with the spheric corona gives the predictions of 
both the disk inclination angle and the absorption that are consistent with the direct 
measurements, and only slightly under-predicts the black hole mass of LMC X-3. 
The model explains the spectral state evolution of LMC X-1 within the BeppoSAX 
observation as a change in the accretion rate, which leads to an increase in both 
the inner disk temperature and the Comptonization opacity. On the other hand, the 
CMCD model with the slab-like corona is more problematic in the test and is thus 
not recommended. 

Key words: black hole physics — stars : individual (LMC X-1, LMC X-3) — X-rays: 
stars 



1 INTRODUCTION 

Accreting black hole X-ray binaries (BH XBs) typically stay 
in on e of the two distinct spectral states jTanaka fc LewinTI 
[^9^. In the so-called high/soft state, such a binary has a 
soft X-ray spectrum and a relatively high X-ray luminosity, 
which is believed to be dominated by the emission directly 
from the accretion disk around the black hole (BH). Under 
the standard geometrically thin and optically thick accre- 
tion disk approximation (see Pringle 1981 and references 
therein), the X-ray spectrum is an integration of the as- 
sumed blackbody-like emission over the disk with a temper- 
ature that decreases with increasing radius jMitsuda et al. I 
Il984l) and is thus called the multi-color disk (MCD) . In the 
opposite low/hard state, the spectrum is relatively flat and 
can often be approximated by a power law (PL) , which may 
extend up to several hundred keV, whereas the luminosity 
is typically low. The flat spectral shape is usually attributed 
to the Comptonization of soft disk photons by hot electrons 
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in a surrounding corona (e.g., Zdziarski 2000). Therefore the 
spectrum of a BHXB system can usually be well fltted by 
a two-component model with a blackbody-like (soft) com- 
ponent plus a PL-like (hard) component; a disk reflection 
component and a broad iron line component are sometimes 
also visible in the spectra of some systems (e.g., Cygnus X- 
1, Di Salvo et al. 2001; GX 339-4, Ueda, Ebisawa, & Done 
1994). 

Of course, there are various other effects that one needs 
to consider. For the soft component, because of the high 
temperature (~ 1 keV) in the inner region of the accretion 
disk, disk emission is expected to be slightly Comptonized by 
the free electrons in this region and the local emergent spec- 
trum could be approximated as a diluted blackbody spec- 
trum rather than a blackbody one. Furthermore, the emis- 
sion from the inner accretion disk is subject to the strong 
gravitational field in the vicinity of the BH, and should also 
be modified by the extreme Doppler motion of the accre- 
tion disk. Detailed calculations suggest that the simple MCD 
model (diskbb in XSPEC) could still be used to describe the 
distorted disk emission, but needs to be corrected by several 
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factors in order to infer the consistent physical quantities: 
a factor f color accounting for the color temperature correc- 
tion of the Comptonization in the inner disk region (e.g, 
Shimura & Takahara 1995; Merloni, Fabian, & Ross 2000); 
a factor r) indicating the difference between the apparent and 
intrinsic radii of the peak temperature, and the integrated 
luminosity difference between the MCD model and a more 
realistic accretion disk with torque-free boundary condition 
(e.g., Kubota et al. 1998; Gierliiiski et al. 1999); factors foR 
and goR accounting for gravitational redshift and Doppler 
shift and the integrated disk flux change due to the general 
relativistic effects (e.g., Cunningham 1975; also see Zhang, 
Cui, & Chen 1997 for more discussions). 

For the hard component, the phenomenological PL 
model is often used. However, the extrapolation of the PL 
model to the lower energy in the MCD-I-PL model neglects 
the seed photon curvature that should be reflected in the 
Comptonized spectrum (e.g., Shrader & Titarchuk 1998; 
Done, Zycki, & Smith 2002). In a typical global fit of the 
MCD-f PL model to a BHXB spectrum, this unphysical ex- 
trapolation usually leads to an artificial increase in the soft 
X-ray-absorbing column density estimates. As a result, the 
absorption-corrected source fiux could be over-estimated. 
This problem is not very acute for an observation with poor 
counting statistics and/or for a low temperature accretion 
disk with the MCD emission peak located outside the ob- 
served energy range (e.g., accreting intermediate-mass BH 
candidates; Wang et al. 2004, Paper II). But for a stel- 
lar mass BH with an accretion disk temperature of ~ 1 
keV, the problem is significant. In fact, it has been shown 
that the MCD+PL model fails to give an acceptable fit to 
X-ray spectra of LMC X-3, when the X-ray absorption is 
tightly constrained by an independent measure ment from X- 
ray a bsorption edges in dispersed X-ray spectra jPaee et al. I 

Instead of using the PL model, several groups have 
developed various Comptonization models by taking care 
of the radiative transfer process either numerically or an- 
alytically (e.g., compbb, Nishimura, Mitsuda, & Itoh 1986; 
comptt, Titarchuk 1994; Hua & Titarchuk 1995; thcomp, 
Zycki, Done, & Smith 1999; eqpair, Coppi 1999 and refer- 
ences therein; compps, Poutanen & Svensson 1996; etc.). 
The seed photon spectrum is assumed to be a single- 
temperature blackbody in compbb or its Wien approxima- 
tion in comptt, which clearly deviates from a multi-color 
black-body disk. This deviation could be very significant, es- 
pecially for a disk with a high inner temperature, typical for 
a stellar mass BH; the black-body would peak well within an 
X-ray band. Whereas the thcomp is a thermal Comptoniza- 
tion model, the eqpair model is a hybrid model of ther- 
mal and non-thermal Comptonization (see also Gierlihski et 
al. 1999); both models can adopt a proper disk spectrum, 
but the latter is more advanced which takes care of nearly 
all the important physical processes in a disk-corona sys- 
tem including Compton scattering, pair production and an- 
nihilation, bremsstrahlung, and synchrotron radiation, etc. 
iCoppi . .1992.) . All these four models can only produce a 
direction-averaged spectrum. Like eqpair model, the compps 
model also contains most physical processes in the disk- 
corona system and allows the use of a MCD spectrum, but it 
can treat geometry more accurately and generate an angu- 



lar dependent spectrum in a spherical or a slab-like corona 
system. 

It is worth noting that in the disk-corona scenario, the 
Comptonization generates the PL-like component at the ex- 
pense of the MCD fiux, therefore the MCD and the PL- 
like components are physically related. Some Comptoniza- 
tion models mentioned above (e.g., eqpair and compps) 
have taken care of this relationship whereas some have not 
(e.g., thcomp). One then should be cautious of interpret- 
ing the physical meaning of the MCD parameters when ap- 
plying the Comptonization models. The MCD model here 
describes only those un- Comptonized (un-scattered or es- 
caped) disk photons rather than the actually original disk 
emission. In particular, this relationship does not exist in 
MCD+PL model, therefore the original disk fiux could be 
under-estimated from the fitted MCD parameters, so are 
other related inferred parameters such as inner disk ra- 
dius. The under-estimate, for example, may be responsi- 
ble for the apparent change in the inner disk radius with 
the transition from one state to another, as has been 
claimed for several sources such as XTE J1550-564, GRO 
J1655-40 ( Sobczak et al. I Il999allb>l . and XTE J2012+381 
JCampana et al. II2OO2I) 

All the models mentioned above have been widely used 
to model the spectra of Galactic BHXBs, neutron star X- 
ray binaries, as well as the extragalactic sources including 
ultraluminous X-ray so urces (ULXs) and AGNs. Ho wever, 
except for the works bv lGierliilski et al. I (liggSllioO j) . there 
has been little rigorous test of the models against direct 
measurements (e.g., neutral absorption column density, BH 
mass, system inclination, etc.), which are available for sev- 
eral well-known systems such as GRO J1655-40, LMC X-1, 
and LMC X-3. 

In Yao et al. (2005; Paper I), we have presented a 
Monte-Carlo method in simulating Comptonized multi-color 
disk (CMCD) spectra. The simulations used the MCD as the 
source of seed photons and self-consistently accounted for 
the radiation transfer in the Comptonization in a spherical 
or slab-like thermal plasma. We have applied this CMCD 
model, implemented as a table model in XSPEC, to a stel- 
lar mass BH candidate XTE J2012-I-381 in our Galaxy. This 
application shows that the inner disk radius is not required 
to change when the source transits from the soft state to 
the hard state, in contrast to the conclusion reached from 
the fits with the MCD-I-PL model dCampana et al. 11200 jl . 
For a spherical corona, the toy model contains the following 
parameters: the inner disk temperature (Tin), the system 
inclination angle (9), the effective thermal electron temper- 
ature (Tc), optical depth (r), and radius (Re) of the corona 
as well as the normalization defined as 

where D is the source distance and Rin is the apparent in- 
ner disk radius. For a slab-like corona, another extreme of 
the geometry, we assume that the corona covers the whole 
accretion disk. We find that the final emerged spectra are 
insensitive to the different vertical scales, so Rc does not 
appear as a parameter in this geometry. Note that Kcmcd 
differs from the normalization for the MCD model. 
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More detailed discussion on the CMCD model can be found 
in Paper II, in which we have applied the same model for 
a spherical geometry to six ULXs observed with XMM- 
Newton. The fitted Tin (~ 0.05-0.3 keV) of these sources 
are distinctly different from the values (~ I keV) obtained 
for known stellar-mass BHs, as presented in this paper and in 
Paper I. Indeed, the inferred BH masses {Mbh) of the ULXs 
are ~ 10^ M©, consistent with the intermediate-mass BH in- 
terpretation of these sources. We have also shown that the 
MCD-I-PL model gives an equivalent spectral description of 
the ULXs, although the CMCD model provides unique con- 
strains on the corona properties and on the disk inclination 
angles, as well as on the BH masses. Because of the lower 
disk temperatures, compared to those of the stellar mass 
BH systems, the nonphysical effects of the MCD-I-PL model 
are typically not significant in the observable photon energy 
range of the intermediate-mass BH candidates. 

In the present work, we conduct a critical test of the 
CMCD and MCD-I-PL models by comparing parameters 
{Mbh, &, and the equivalent neutral hydrogen absorption 
Nh) inferred from the X-ray spectra of LMC X-1 and X- 
3 with the more direct measurements based on optical and 
dispersed X-ray spectra; we also compare the results from 
the different corona geometrical configurations of the CMCD 
model. We first briefly describe these measurements and the 
X-ray observations in §2, and then present the spectral fit- 
ting results in §3. We describe the specific comparisons in 
§4 and present the discussion and our conclusions in §5. 



2 DESCRIPTION OF THE SOURCES AND 
OBSERVATIONS 

We select LMC X-1 and X-3 for this study chiefly be- 
cause of their location in our nearest neighboring galaxy, 
the Large Magellanic Cloud (LMC; D = 50 kpc is adopted 
throughout the work). Both the well-determined distance 
and the relatively low foreground soft X-ray absorption are 
essential to our test. These two sources are also among the 
three well-known persistent BHXBs and are usually found 
in the high/soft state. The low /h ard state was occasionally 
repor ted for LMC X-3 llBovd fc Smale .2000. : Homan et al. I 
I200CI) . but never for LMC X-1. The remaining known persis- 
tent BHXB, Cygnus X-1, is in our Galaxy and stays mostly 
in the low/hard state (e.g., Pottschmidt et al. 2003). The X- 
ray spectra of this source also show a strong disk reflec tion 
component JCilfanov et al. Ill999t [Prontera et al. 11200 — 
a complication that is not included in the CMCD models. 
Table summarizes the key parameters of LMC X-1 and 
X-3, which are used for the comparison with our spectrally 
inferred values (§4). 

Both LMC X-1 and X-3 were observed with Chan- 
dra (e.g., Cui et al. 2002) and XMM-Newton (e.g.. Page 
et al. 2003). The dispersed X-ray spectra of LMC X-3 have 
been used to measure the X-ray absorption edges (mainly 
for Oxygen) , which tightly constrains the the ab sorbing mat- 
ter co lumn density Nh along the line of sight llPaee et al. I 
|20o3). The absorption towards LMC X-1 is, however, sub- 
stantially higher. As a result, the photon flux at the Oxygen 



edge is too low to allow for a useful constraint on Nh based 
on the existing data. 

We here utilize the data from the BeppoSAX observa- 
tions, which were carried ou t on 1997 October 5 for LMC 
X-1 and October 11 for X-3 iTreves et al. II2OO0I) . The data 
do not have pile-up problems, which could be present for 
X-ray CCD imaging observations of bright sources. Four 
types of narrow-field instruments (NFIs) were on board: 
Low Energy Concentrator System (LECS), Medium Energy 
Concentrator Systems (MECS), High Pressure Gas Scintil- 
lation Proportional c ounter (HPGSPC), and Phoswich De- 
tector System (PDS) jBoella et al. Ill997^ . The exposure for 
LMC X-1 and X-3 were about 15 ks each for the LECS 
and about 40 ks each for the MECS. These two instruments 
were sensitive to X-rays in the energy ranges of 0.1-10 and 
1.3-10 keV respectively. Data from the HPGSPC and the 
PDS, which were sensitive to photons in 4-120 keV and 15- 
300 keV ranges respectively, were not included because of 
poor counting statistics, and also because of possible source 
contaminati on from PSR 0540-69 (which is 25' away from 
LMC X-1) JSeward et al. Ill984l : iHaardt et al. ll200lD . We 
extracted the spectra from a radius of 8' and 8.4' around 
each source from the LECS and the MECS observations 
and used the energy ranges of 0.2-4 keV and 1.8-10 keV for 
these two instruments in this study. The background con- 
tributions to the spectra are small and are estimated from 
a blank field. The spectra from the LECS, the MECS2, and 
the MECS3 were jointly fitted for each source, using the 
software package XSPECn.2.0bs. 



3 RESULTS 

Wc summarize the spectral fitting results and the inferred 
source fiuxes in Table |21 The quoted uncertainty ranges of 
the parameters are all at 90% confidence level. Fig. shows 
the spectral fits with the CMCD models. The systematic de- 
viation of the data from the model at low energies (< 1 keV) 
might bo duo to poor ca libration of the instrument spectral 
response ( Martin et al. ]il996). Fig. H illustrates the effects 
of the Comptonization in the spherical corona systems. 

Both CMCD and MCD-I-PL models give acceptable 
fits. The model parameters are all well constrained except 
for the CMCD parameter 6, for which only the upper or 
lower limit is constrained. The MCD-I-PL model parame- 
ters we obtained her e are consistent with those reported by 
iHaardt et aTl J200ll) . The fitted Nh values from MCD-I-PL 
are systematically higher than those from CMCDs. The 
same is true for the inferred absorption-corrected fluxes, es- 
pecially for LMC X-1 which is a factor of ~ 4 higher from 
MCD-hPL than from CMCD (Table HJ. 

For LMC X-3, Ti„ from MCD-^PL is sUghtly higher 
than those from CMCDs, whereas the value of Kmcd /cos{9) 
is consistent with those of Kcmcd- The best flt parameters 
in the two different geometric CMCDs are nearly identical. 
The small value of Tc indicates that only a small portion of 
disk photons have been up-scattered to high energies. 

For LMC X-1, except for a consistent Nh value, the 
fitted parameters are significantly different between the two 
different geometric CMCDs (Table 0. We wiU see in §4, 
the results from the slab-like configuration are inconsistent 
with the independent measurements. Tin values from the 
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Table 1. Comparison of parameter measurements. 





e 






T 






n 


(Me) 


(lO^fcm^^) 


(keV) 


References 






LMC X-1 






Indep. est. 




4^M^12.5 






1, 2 


CMCD: sphere 


<43 


4.0(3.8-4.5) 


54(52-56) 


0.93(0.91-0.94) 




CMCD: slab 


>64 


15.2(10.4-19.9) 


53(51-54) 


0.80(0.79-0.81) 




MCD+PL 


N/A 


3.0(2.9-3.1) 


84(79-89) 


0.93(0.92-0.95) 








LMC X-3 






Indep. est. 


e <7o 


M> 5.8 ±0.8 


3.8(3.1-4.6) 




3,4,5,6 


CMCD: sphere 


<69 


4.19(4.17-4.21) 


4.5(4.2-4.7) 


0.98(0.97-0.99) 




CMCD: slab 


<61 


3.73(3.71-3.76) 


4.4(4.2-4.6) 


0.96(0.95-0.97) 




MCD+PL 


N/A 


4.2(4.1-4.3) 


7.0(6.0-8.0) 


1.02(1.01-1.03) 





T he BH mass VLbh is e stim ated by assumin g zero spi n of the BH. S ee tex t f or details. Referen ces: 
1 'Hutchi ngs et al I ^1983^: ^ iG icrlinski et al. I <2001^ : ^ ICowlev et al. I IWHct) : IPaczvnski I IWHct) : ^ 
ISoria et al. I ^2001^ : "^ iPage et al. . I.200a'l . 



Table 2. Spectral fit results 




1 10 1 10 

Channel energy (keV) 

Figure 1. Model fits to the BcppoSAX spectra of LMC X-1 
{panel a) and LMC X— 3 {panel c), and the corresponding residuals 
in term of sigmas [panels b and d). The solid iine in panels a 
and c show the fit of CMCD model. The fit goodnesses from the 
different geometrical models (a sphere vs. a slab) arc nearly the 
same. 
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Figure 2. The effects of the Comptonization in LMC X— 1 (a) 
and in LMC X-3 (b). Solid line: spherical CMCD model with 
the best fit parameters; dotted line: CMCD model with the same 
parameters but t=0, which is equivalent to the MCD model. 



parameters 


LMC X-1 


LMC X-3 




CMCD: sphere 




Nh (1021 cm2) 


5.4(5.2-5.6) 


0.45(0.42-0.47) 


(keV) 


0.93(0.91-0.94) 


0.98(0.97-0.99) 


(keV) 


19(15-23) 


20(19-23) 


Rc (Rg) 


11(9-19) 


10.0(9.6-16.6) 


T 


1.0(0.8-1.3) 


0.10(0.09-0.13) 


e (deg) 


28(**-43) 


59(**-69) 


-f^CMCD 


57(51-73) 


44.0(43.5-44.3) 


xVdof 


641/627 


760/649 


/0.2-10 


8.6 


5.7 




CMCD: slab 




Nh (10^1 cm2) 


5.3(5.1-5.4) 


0.44(0.42-0.46) 


Ti„ (keV) 


0.80(0.79-0.81) 


0.96(0.95-0.97) 


(keV) 


10.0(9.0-10.6) 


22(19-24) 


r 


0.45(0.41-0.49) 


0.09(0.08-0.12) 


e (deg) 


75(64-**) 


50(**-61) 


-f^CMCD 


489(230-840) 


42.7(42.1- 43.3) 


xVdof 


638/628 


759/650 


/0.2-10 


8.5 


5.7 




MCD+PL 




Nh (10^1 cm2) 


8.4(7.9-8.9) 


0.7(0.6-0.8) 


Tin (keV) 


0.93(0.92-0.95) 


1.02(1.01-1.03) 


-f^MCD 


28(26-30) 


23(22-24) 


r 


3.5(3.4-3.6) 


2.6(2.5-2.8) 


KpL (10-1) 


2.3(1.9-2.7) 


0.19(0.14-0.23) 


xVdof 


624/629 


747/651 


/0.2-10 


31 


6.4 



The uncertainty ranges are given in parenthesis at the 90% con- 
fidence level; asterisks indicate that the limit is not constrained. 
Rg = GM/c'^. The /0.2— 10 is the absorption-corrected flux in the 
energy range 0.2-10 keV and in unit of 10" i*^ erg cm" ^ s^i. 

spherical CMCD and the MCD+PL are consistent with each 
other, but Kmcd/cos{9) value, assuming the best-fit CMCD 
6, is ~ 1.5 smaller than Kcmcd- 

Because a significant change in the spectral shape oc- 
curred during the observation of LMC X-1, we have fur- 
ther split the exposure into two parts, the first 30 ks and 
the re maining time, in the same way as in iHaardt et al. I 
i200ir) . To tighten the constraints on spectral parameters, 
we jointly fit Nh and 9, which should be the same in the 
two parts of the observation (Table EJ. From the early part 
to the later part, according to the spherical CMCD model. 
Tin increased by ~ 7% and the corona became a factor of 
~ 2 larger and opaque with r increased by a factor of ~ 4. 
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Table 3. Spectral variation of LMC X-1 



parameters 


part 1 


part 2 




CMCD: sphere 




Nh (lO^i cm2) 


5.3(5.2-5.5) 


= part 1 


(keV) 


0.91(0.90-0.92) 


0.96(0.93-0.99) 


(keV) 


13(9-18) 


10(8-12) 


Rc (Rff) 


10(9-16) 


25(18-33) 


T 


0.55(0.5-0.7) 


2.0(1.8-2.4) 


6 (deg) 


23(**-45) 


= part 1 


KcMCD 


61(55-72) 


57(46-77) 


X^/dof 


1211/1192 


/0.2-10 


8.8 


8.3 




CMCD: slab 




Nh (10^1 cm^) 


5.3(5.2-5.4) 


= part 1 


(keV) 


0.81(0.80-0.82) 


0.74(0.73-0.76) 


Tc (keV) 


20(19-23) 


5(**-7) 


T 


0.11(0.10-0.14) 


1.9(1.7-2.1) 


(deg) 


75(64-**) 


= part 1 


KcMCD 


348(224-495) 


1223(730-1855) 


xVdof 


1209/1194 


/0.2-10 


8.8 


8.1 




MCD+PL 




Af_F/ (10^1 cm^) 


8.5(8.0-9.0) 


= part 1 


(keV) 


0.90(0.89-0.91) 


1.00(0.98-1.03) 


Km CD 


52(49-55) 


23(21-26) 


r 


3.6(3.5-3.8) 


3.4(3.3-3.6) 


^PL (10-1) 


3.0(2.4-3.7) 


3.4(2.8-3.9) 



x'^/dof 1180/1195 

/0.2-10 48 42 

Please refer to Tabic |2l 

The source flux in the 0.2-10 keV band decreased slightly, 
although the normalization remained essentially the same. 
Because Ti„ oc (M /Mbh)^^'^ , where Mbh and M are the 
BH mass and the accretion mass rate, a rising Ti„ during the 
observation of LMC X-1 was then caused by an increasing 
accretion rate. Apparently, this change led to the thickening 
of the corona. The slight decline of the flux is likely due to a 
combination of the energy loss to the corona and to the scat- 
tering of photons to energies greater than > 10 keV. Fig. |3] 
demonstrates the differences of the Comptonization effects 
in the two parts of the observation. However, it is hard to 
physically understand the fitted parameters of the slab-like 
CMCD: From early part to the later part, Ti„ and Tc need to 
decrease by ~ 16% and a factor of 4, respectively, but in the 
mean time both r and Kcmcd (hence Rin) have to increase 
dramatically. We believe that the slab-like CMCD model is 
not suitable for describing such a spectral state of a BHXB 
system when the hard component contributes significantly 
to the total flux. 

Now let us check the results from the M CD-I-PL model, 
which are fairly consistent with those of iHaardt et al. I 
i200ll) . The MCD normalizations are significantly differ- 
ent between the early and later pa rts of the observatio n 
(Table 1^, confirming the analysis bv lHaardt et al. I J200ll) . 
which is based on the same data and same model. They 
concluded that as Tin increased (by ~ 9%), Rin decreased 
(~ 38%) from the early to later parts of the observation. 
But, as was discussed in §1, this apparent change in Rin is 
most likely due to the lack of the accounting for the radia- 
tive transfer between the two components of the MCD -l-PL 
model and is therefore not physical. 

Furthermore, Fig.^Jshows clearly that for both sources. 



1 10 
Channel energy (kcV) 

Figure 3. The spectral variation of LMC X— 1 from the first 30 
ks (part 1, solid line) to the remaining time (part 2, dashed line). 
The models are plotted with the best fitting parameters (TableUJ. 
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Figure 4. The comparison of the spherical CMCD model and 
MCD-I-PL model with the best fit parameters (Table |5) for (a) 
LMC X-1 and (b) LMC X-3. 

the PL component surpasses the MCD component in con- 
tributing to the spectra at low energies ( < 1 keV for 
LMC X-1 and < 0.4 keV for LMC X-3). This nonphysical 
straight extension of the PL component to the low energy 
parts of the spectra is the main cause for the required high 
values of Nh in the MCD-I-PL fits, compared to those in the 
CMCD fits (Tables H and EJ. 



4 COMPARISONS WITH INDEPENDENT 
MEASUREMENTS 

In addition to the above self-consistency check of the X-ray 
spectral models, we compare the inferred parameter values 
of the BHXBs with the independent measurements to fur- 
ther the test. The key results of this comparison are included 
in Tableland are discussed in the following: 
• X-ray Absorption Nh may be obtained more directly 
and accurately via the spectroscopy of neutral absorption 
edges of Oxygen and Neon (e.g., Schulz et al. 2002; Page et 
al. 2003). This method is nearly independent of the over- 
all continuous spectral shape and any other source prop- 
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erties. Usin g the dispersed spe ctrum of an XMM-Newton 
observation. jPaee et aTl J2003h inferred Nh = O.SStoioy x 
10'^^ cm~^ for LMC X-3, as s uming the intersteUar medium 
abundance of IWilms et al. I ll200Cf) . We have also adopted 
this assumption by using the absorption model TBabs in 
XSPEC. The Nh values from the CMCD models and from 
the X-ray absorption edge measurement agree with each 
other within the quoted error bars, whereas the value from 
the MCD+PL model is significantly higher (Table Q. 

• Disk Inclination Angle The upper limits on 9 from 
the spherical CMCD model are consistent with the values 
obtained from the optical observations (Table 0. There is 
a simple reason why 9 can be constrained in the spherical 
CMCD model: the Comptonized flux is nearly isotropic and 
barely affected by the disk inclination, whereas the observed 
strength of the soft disk component is proportional to cos(6'). 
This mostly geometric effect is strong when 9 is large; e.g., 
no radiation come directly from an edge-on disk. Therefore, 
the upper limit can be constrained reasonably well, which is 
especially important for estimating the BH mass (cf. Eqs.0 
and|21l. In the slab- like CMCD model, since both of the ob- 
served soft and hard components strongly depend upon the 
system inclination angle (e.g., Sunyaev & Titarchuk 1985; 
see also the Fig. 2 in Paper I), 9 in principle can be well 
constrained in this model. For LMC X-3, the PL-like com- 
ponent only contributes a small portion of the total fiux, 
and the constraint of 9 mainly from the soft component, 
as in a spherical CMCD model. For LMC X-1, the PL-hke 
component contributes significantly to the total flux. The 
slab-like corona configuration becomes problematic, giving 
the inconsistent 9 constraints. 

• Black Hole Mass The mass of each putative BH may 
be estimated as M = f Rin/Ga, where a = 6 or 1 for a 
non-spin or extreme spin BH), fRin is the inner disk radius, 
assumed to be the same as the radius of the last marginally 
stable orbit around the BH, and the factor / [depending 
on the system inclination, 0.94 for LMC X-1 and 1.12 for 
LMC X-3; we adopt a new valu e of 77 = 0.41 (see §1) de- 
rived from lKubota et al. J l| l998^. compared to the old val- 
ues used 77 = 0.7 used in HhariE et al. (1993) and in Paper 
II] includes various corrections related to the spectral hard- 
ening, special and general relativity effects (e.g., Cunning- 
ham 1975; Zhang, Cui, & Chen 1997; Gierlihski et al. 1999, 
2001). Assuming no spin for the BHs, from the spherical 
CMCD results, we estimate the BH masses as 4.0(3.8-4.5), 
4.19(4.17-4.21) Mq for LMC X-1 and LMC X-3, respec- 
tively. For LMC X-1, the value is consistent with the result 
from the optical study, whereas for LMC X-3, the derived 
Mbh is slightly smaller (Table which may suggest that 
the LMC X-3 is a mild spin system. 



5 SUMMARY 

In this work, we have applied the MCD-I-PL model as well 
as our recently constructed CMCD model to BHXB systems 
LMC X-1 and LMC X-3, confronting the fitted parameters 
with directly measured values. We have also tested two dif- 
ferent corona geometric configurations. The spherical con- 
figuration passes almost all the tests: the effective hydrogen 
column densities, disk inclinations, and the BH masses of 
the LMC X-1 and LMC X-3. This consistency suggests that 



the CMCD model with a spherical corona provides a reason- 
ably good spectral characterization of BHXBs. The model 
offers useful insights into physical properties of the Comp- 
tonization coronae and their relationship to the accretion 
process. In contrast, the slab-like CMCD model is problem- 
atic in describing the spectrum of LMC X-1, in which the 
PL-like component contributes significantly. Similarly, the 
MCD-I-PL model, though generally providing a good fit to 
the spectra of BHXBs, could give misleading parameter val- 
ues. Although the tests conducted in this work are still very 
limited, they have demonstrated the potential in discrimi- 
nating among various models. 
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